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Abstract 

We present an electron projection imaging method to study the ultrafast evolution of photo- 



electron density distribution and transient fields near the surface. The dynamical profile of the 
photoelectrons from graphite reveals an origin of a thermionic emission, followed by an adiabatic 
process leading to electron acceleration and cooling before a freely expanding cloud is established. 
The hot electron emission is found to couple with a surface charge dipole layer formation, with a 
sheet density several orders of magnitude higher than that of the vacuum emitted cloud. 
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Understanding the mechanisms of vacuum space charge (VSC) 

HQS 

emission and 

surface charge formation is central to the development of pulsed laser driven electron tech- 
nologies such as time-resolved pho toemission AS , |4|, scanning probe microscopy 



electron diffraction 



b, y, y, q, 
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12j , and micro scop 



charge has been studied both theoretically 
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13l . |14| . While the vacuum space 
18j] and experimentally 



the role of subsurface electron dynamics within the materials and the nature of the early 
development of VSC remain topics of high interest. Recently, the transient vacuum electric 
field established by the generation of VSC has been investigated, based on the influence of 
the field on a pulsed electron beam 20|, |21j. Here, a method to directly image the spa- 
tiotemporal evolution of the photo-emitted electron bunch generated over a femtosecond 
laser excited surface is presented, based on an electron projection geometry. The method 
possesses sufficient sensitivity to image electron bunches as small as 10 10 e/cm 3 and permits 
quantitative measurement of the instantaneous electron bunch density distribution and its 
translational and expansion velocities in the picosecond and micrometer regime. Contrary 
to the space-charge heating, we observe an adiabatic cooling during the initial expansion 
of the thermally emitted electrons, causing a nearly 80% drop of the internal temperature 
while accelerating a high CoM velocity. In conjunction with the diffractive voltammetry 



121 . |22j , this offers a way to study simultaneously both subsurface charge dynamics as well 
as vacuum emitted space charge effects self-consistently. 

The ultrafast electron imaging is conducted in a pump-probe experimental arrangement, 
in which the femtosecond laser (800 nm, 45 fs) is used to generate the photoemission from 
graphite (HOPG, ZYA grade, SPI Supplies), which is subsequently probed by the delayed 
surface scattered electrons. The incident laser beam is displaced from the source of electron 
scattering (P) by a distance xq (typically a few mm) to establish a projection imaging ge- 
ometry (Fig. [1]). The presence of VSC diminishes the forward-scattered electrons generated 
at P, thus casting a shadow on the screen, which is at a distance L = 150 mm away from 
P. The magnification from this projection imaging is M m L/xq ~ 100. By adjusting the 
arrival time of the probe electrons relative to that of the laser, shadow images of the evolving 
electron bunch can be obtained, as shown in Fig. 2. The cross-sectional line scans obtained 
from these shadow images reveal an accumulation peak near the shadow edge (caused by 
surface scattering), which decays sharply into the vacuum, followed by a depletion profile 
characteristic of a gaussian evolving in space over time. To extract the spatiotemporal dy- 
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FIG. 1: (Color online) Schematics of the electron point-projection imaging technique. Typical 
values are xq ~ 1 - 4 mm, L ~ 150 mm. The dashed arrows illustrate the effect of the photo-emitted 
electron bunch on the Bragg beams. 

namics of the electron bunch, we fit the line scans to the following analytical form, which 
takes into account the effect of the projection geometry: 

/ (dxp-Lzp) 2 \ 

F(d) = A exp(-d/A ) + sf-exp - ffl±£f (1) 

where d is the distance on screen measured from the shadow-edge. The charge distribution 
p(z, t) is modeled as gaussians in x, y and z directions with a sheet electron density So in 
the xy-plane, 1/e half-widths A^, A y , A z and a CoM position (x Q , y , z ). A x « 500 jj,m 
is determined by examining the transverse size of the electron bunch, which remains nearly 
unchanged and corresponds well with the width of laser illumination on the surface. The 
evolution of the bunch in the z-direction is shown in Figs. [3] (a), (b) and summarized in 
Table I for ffuences F =7, 23 and 56 mJ/cm 2 . The results indicate the electron cloud CoM 
to follow a projectile-like trajectory (caused by the attractive force from image charges, Fig. 
El^a)), with a translational CoM velocity Vq = dz /dt, while undergoing free-expansion as 
indicated by the nearly constant expansion speed v e = dA z /dt in Fig. E](b). 

From the gaussian shape of the velocity distribution and the lack of a power-law enhance- 
ment of S over increasing F, we believe a thermionic emission scenario is best to describe 
the dynamical parameters observed here. For a steady-state thermal emission, the initial 
translational speed Vi and the expansion speed v e of the bunch is related to the electronic 




FIG. 2: (Color online) (a) Snapshot shadow images of photoemitted electron bunch obtained with 
xq = 4.3 mm. The cross-sectional line profiles (black curves) are extracted along the dashed line 
in the second panel, (b) Fitting curves (lines) of cross-sectional line-profile data (circles) using a 
gaussian model considering the projection geometry. 

temperature (T e ) at the surface via the relations: Vi = <J kBT e /2iim e , and v e = J kBT e /m e . 
The observation of a vq significantly higher than v e in Table [I] suggests an electron accel- 
eration and cooling during the initial adiabatic expansion of the cloud. This scenario is 
supported by the observation of an apparent shift of the zero-of-time extracted from the 
linear extrapolation of cloud expansion to intercept with the time axis (Fig. 3(b)), indicat- 
ing a rapid decrease of v e within 10 /im (resolution-limited) of the expansion. This cooling 
process converts part of the electronic enthalpy (Ah = 5/2/csAT) to the kinetic energy 
of the bunch, and under such an adiabatic model we can deduce the initial electron tem- 
perature (Tj) using l/2m e (Vg — vf) = Ah. The Tj deduced from our experiment shows a 
saturation at high fluence, and the highest electronic temperature obtained is just above 
the laser pulse energy 1.55 eV (18,000 K), further confirming a thermally limited emission. 
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FIG. 3: (Color online) Photoemission from graphite surface: Dynamics along the z-direction, 
when xq = 2.4 mm. (a) Time-evolution of the electron bunch's center-of-mass (CoM) and (b) 
1/e- width, (c) Spatiotemporal evolution of the Bragg peaks and the electron bunch at F = 56 
mJ/cm 2 , showing the bunch CoM (black dots) and its spread (1/e- width) relative to the CoM 
(gray circles). 

During the same period, by assuming an initial cloud size of 30 nm (comparable to the laser 
penetration depth) and Eo = 2 x 10 8 e/cm 2 , we estimate the possible heating effect arising 
from space-charge driven acceleration of the electrons at the expense of Coulomb self-energy 
to be at most 0.2 eV. Thus Coulomb explosion cannot provide the observed high initial CoM 
velocity. 

In addition to the imaging, we also analyze the trajectory of the diffracted electron 
beams (0,0,6), (0,0,8), (0,0,10) in the central streak region of the diffraction pattern. The 
instantaneous vacuum electric field E(z,t) established by the two opposing fields associated 
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TABLE I: Results of fitting experimental data to the projected gaussian model. 
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with the VSC p(z, t) and its mirror- charges on the surface can be described by : 

E(z, t) = ^- (a(z)E - 2 Jp(z', t)dz'\ (2) 



where a(z) z/yz 2 + A 2 is the proximity effect factor [201] . describing the reduction in 
local electric field from mirror charges with a finite-sized slab, which is [l-a(z)]eSo/(2eo)]. 
The Bragg beams located at a height z c above the surface interact with VSC and exhibit 
an inversion, as shown in Fig. [3](c) when the CoM of the electron bunch reaches a height 
z o — z c~l^zi such that the local electric field at z c , E(z c ) = 0. Carrying out the integration 
in Eqn. El we obtain the condition for Bragg beam inversion as 7 = 1.6, 1.3, 1.1 for zq = 
50 /im, 100 /im and 150 \im corresponding to the (0,0,6), (0, 0, 8) and (0, 0, 10) diffracted 
beams respectively. Inspecting the trajectories A z (t) and zo(t) in Fig. 3(c) at the point 
of Bragg beam inversion, we measure 7 = 1.0 ± 0.4, 1.2 ± 0.2 and 1.2 ± 0.2, which is 
in good agreement with the predicted values above. The deviation in the case of (0,0,6) 
beam is believed to be caused by the recoil effect, which plays a more significant role here. 
Given this confirmation, we can reliably obtain the sheet electron density So for the three 
fluences studied here (Table I), and compare them with the initial electron temperature Tj. 
S scales linearly with Tj to a satisfactory degree, agreeing with the thermionic emission 
model proposed by Downer and coworkers 231 ] . The electron bunch dynamics measured here 
reveals that the majority of the electrons will return to the surface within 1 ns, see t return in 
Table I, with only a very small portion (< 10 -4 ) of the cloud remaining in vacuum beyond 
3 ns. 

We also observe a dramatic increase in the electron refraction {56b) [22] simply shifting 
the laser to P without changing the diffraction geometry, as shown in Fig. H^a). Simulation 
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FIG. 4: Comparision of TSV and VSC effects in UEC investigations of HOPG. (a) Shift of 
the (0,0,6) Bragg peak in case of perfect pump-probe alignment (xq = 0) and projection imaging 
geometry with (xo = 2.4 mm), (b) Comparison of the TSV measured by the probing electron beam 
submerged beneath the HOPG surface in the UEC study (Ref. 24|) along with the estimated effect 
of VSC in that study, based on the VSC dynamics extracted here. 

of the VSC induced refraction effect in the case of perfect pump-probe alignment (Fig. |4^b)) 
clearly indicates that VSC alone is not sufficient to explain the observed large and rapid 
10 V rise in the transient surface voltage of fs-laser excited graphite surface 24j. This 
consequently mandates the presence of a surface dipole field that is invisible in vacuum to 
account for the full refraction effect. Assuming a probe-electron depth of ~1 nm, we deduce 
a ps retention of surface charge density in graphite on the order of 10 14 e/cm 2 , which is 
comparable to the typical surface/interface state density. The thermionic emission on a sub- 
to a few picosecond timescale is likely mediated by the image/interlayer states of graphite, 
which has a strong 3D character!^. These states above the vacuum level can efficiently 
transport electrons to vacuum and create the dipole field near the surface. Following the 
electron cooling, the coherent interlayer transport normal to the basal planes of graphite is 
essentially turn off, causing the slow decay of the surface dipole field, which plays a major 
role in influencing the surface structure dynamics 24, Q, [27J]. In conjunction with transient 



surface voltammetry, this vacuum electron imaging method can be applied to investigate a 
variety of hot electron processes involving intense laser interaction with solid surfaces. 
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